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Abstract 

We further develop an extended dynamical mean field approach introduced 
earlier. It goes beyond the standard D = oo dynamical mean field theory 
by incorporating quantum fluctuations associated with intersite (RKKY-like) 
interactions. This is achieved by scaling the intersite interactions to the same 
power in 1/D as that for the kinetic terms. In this approach, a correlated lat- 
tice problem is reduced to a single-impurity Anderson model with additional 
self-consistent bosonic baths. Here, we formulate the approach in terms of 
perturbation expansions. We show that the two-particle vertex functions 
are momentum-dependent, while the single-particle self-energy remains local. 
In spite of this, the approach is conserving. Finally, we also determine the 
form of a momentum-dependent dynamical susceptibility; the resulting ex- 
pression relates it to the corresponding Weiss field, local correlation function 
and (momentum-dependent) intersite coupling. 

PACS numbers: 71.10. Hf, 71.27.+a, 71.28.+d, 74.20.Mn 
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I. INTRODUCTION 



In strongly correlated electron systems, both local and non-local interactions are impor- 
tant in determining the nature of the ground state and low lying excitations. One example 
illustrating this point comes from Kondo systems, such as heavy fermions. Here, the compe- 
tition between the local Kondo interactions and non-local RKKY interactions was recognized 
to play an essential role from very early onEH. The Kondo effect tends to quench local mo- 
ments altogether, while the RKKY interactions promote magnetic ordering. What happens 
when the two processes are about "equally" important is an intriguing question which re- 
mains poorly understood. This question has once again become centrally important, due to 
the emergence of a host of heavy fermion metals lying in the vicinity of a quantum phase 
transition!!. 

The interplay between local and non-local interactions is also essential in Mott-Hubbard 
systemsB. When on-site interactions are strong, their effects can be thought of as determining 
the atomic configurations that lie at low energies. The precise form of the ground state and 
low-lying excitations, on the other hand, have to be determined by the residual intersite 
couplings between these low energy configurations. 

In general, a separation of electron-electron interactions into local and non-local ones is 
not necessarily well-defined. Local interactions, when combined with kinetic terms, can lead 
to effective non-local interactions. After all, both RKKY and super-exchange interactions 
arise in this fashion. 

In theoretical approaches, however, such a separation can often be sharply defined. In 
this paper, we are concerned with the dynamical mean field theory (DMFT)i which is 
formally exact in the limit of infinite dimensions (D = oo)i. The DMFT reduces a correlated 
lattice problem to a self-consistent Anderson impurity model, namely a quantum impurity 
coupled to a self-consistent fermionic bath. The interactions between the impurity degrees 
of freedom reflect the on-site interactions of the lattice problem; in this way local quantum 
fluctuations are retained. The self-consistent fermionic bath of the impurity problem reflects 
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the influence, at the one-particle level, of the rest of the lattice on the selected (i.e., impurity) 
site. All the intersite correlations of the lattice problem, on the other hand, are neglected. 
In this sense, non-local quantum fluctuations are completely lost. 

In earlier works, wefli and independently Kajueter and Kotliari have extended the the 
DMFT such that intersite quantum fluctuations are treated on an equal footing with local 
ones. This extended DMFT reduces a correlated lattice problem into a novel effective 
impurity problem, which corresponds to an Anderson impurity model with additional self- 
consistent bosonic baths. These bosonic baths reflect the influence, at the two-particle level, 
of the rest of the lattice on the impurity site. Through self-consistency, they keep track of 
the intersite quantum fluctuations. In these earlier works, the mean field equations were 
derived using the so-called "cavity" -method. 

The purpose of this paper is to give an alternative formulation of this extended DMFT us- 
ing perturbation methods. We are then also able to explicitly show the conserving nature of 
this approach. Finally, we also derive the expressions for momentum- dependent correlation 
functions, which specify how to calculate the correlation functions from the corresponding 
Weiss fields, local correlation functions and (momentum-dependent) intersite interactions. 

To be specific, we focus on a one band model, 

H = Y^Un^n iX + Ujcl a c ja 

i <ij>,a 

+ o v ij( n i~ < n >)( n j- < n >) + \ J ifii ' 4 (!) 

<ij> <ij> 

The first two terms alone would correspond to the standard Hubbard model for a spin 1/2 
band. The third and forth terms are the intersite density-density (vij) and spin-exchange 
(Jij) interactions. Here n« and Si are the density and spin operators for the c— electrons. 
< ij > labels a pair of nearest-neighbor sites. For simplicity, we limit both the hopping and 
intersite-interaction terms to nearest-neighbor only. Generalizing to the case with longer- 
range hopping and interaction terms is straightforward. 

In the standard large D approach, the single-electron hopping term is taken to be of 
order 1/yD: 
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t Kij> = t/y/D (2) 

The large D limit leads to an effective single-site problem, which describes an impurity 
coupled to a self-consistent Weiss field. There is one Weiss field for each frequency, due 
to the quantum mechanical nature of the hopping term. The coupling of the impurity to 
this frequency-dependent Weiss field can be equivalently described in terms of a coupling 
between the impurity and an effective non-interacting fermionic bath. The bath is fermionic, 
since the Weiss field describes the influence of the rest of sites to the selected impurity site 
at the one-particle level. On the other hand, the intersite interaction terms (both v<ij> and 
J<ij>) are taken to be of order 1/D. With such scaling, only the static Hartree contributions 
survive the large D limhM All the quantum fluctuations associated with these interactions 
are then neglected. 

In the extended DMFT, these intersite interaction terms are also scaled to the order 
1 / \f~D. The large D limit then leads to a different effective impurity problem: the impurity 
is now also coupled to the frequency-dependent Weiss field induced by these intersite interac- 
tions. The two-particle nature of the intersite interactions dictate the bosonic nature of the 
corresponding effective baths in the impurity problem. As a result, the effective single-site 
problem can be thought of as an impurity coupled not only to a self-consistent fermionic 
bath but also to self-consistent bosonic baths. 

The rest of the paper is organized as follows. In Section ||, we introduce the extended 
large D limit and derive the mean field equations using perturbation methods. We then 
derive the expressions for the momentum-dependent correlation functions in Section PH . 
The proof that the approach is conserving is given in Section [TV]. In Section [V], we a) 
generalize the approach to multi-band systems as well as to the case of an ordered state; 
b) specify an approximate procedure to deal with incommensurate spatial fluctuations; and 
finally c) compare our approach with others within the general DMFT framework. 
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II. THE EXTENDED DYNAMICAL MEAN FIELD APPROACH 



In this section, we introduce the extended large D limit and derive the mean field equa- 
tions. 

The intersite hopping term, £<«.,■>, remains of order 1/y/D as given in Eq. (0). The 
nearest-neighbor interactions are now scaled to the same order, 

v<ij> = vj VD 

J <VJ> = J/VD (3) 

In order for the large D limit to be well-defined, we need to subtract the Hartree contribution 
as has already been done in Eq. ([!]). 

We now establish that the self-energy is still local. Consider first the self-energy £<y>. 
The Hartree contributions from both v < ij > and J<ij> vanish. The Fock and high-order 
contributions to the self-energy can be written in a skeleton expansion. As illustrated in 
Fig. H|, any skeleton expansion diagram for the self-energy contains at least an intersite 
interaction path and a fermion propagator from site % to site j. Both are at least of order 
1/\AD. Therefore, S <ij> ~ 0(1/ D). More generally, Ey ~ 0(1/D)^~^^, where ||z— j\\ is the 
least number of steps from site % to site j. This implies that the self-energy is momentum- 
independent: E(q, u) = Ejj(a>). 

Consider now the on-site self-energy Hh(uj). The only real-space self-energy diagrams 
that survive the large D limit have the form illustrated in Fig. [|a), as explained in detail in 
Appendix 0. Here a solid line represents the fermion propagator, Gu. A dashed line denotes 
the intersite interaction (either Vij or Jij). The loop formed by two dashed lines enclosing 
a solid square represents either x^o^) or X^oty) which, as explained in Appendix 0, has 
the following form, 

,ij Xch,iOXch,Oj/ Xch,Oo) 
Xs,0 = ^2 JioJ()j(Xs,ij — Xs,i0Xs,0j/Xs,00) (4) 
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where Xch and Xs are the charge and spin susceptibilities respectively. 

The above implies that can be equivalently calculated in terms of a local problem 
with an action of the following form, 

S MF = f dr Un^n^r) - f dr f dr' \£ d{r)G^ \r - r')c a (r') + 

JO JO Jo a 

+ ■ n(r) : X ^ (t - r') : n(r') : +S(r) ■ Xs J(r - t')S(t')] (5) 

where : n := n— < n >. The skeleton expansion for the self-energy of this local problem, 
S ioc , has the form given in Fig. |]b) where a solid line represents G\ oc and a shaded line 
represents either x7ho or X^o- Since there exists a one-to-one correspondence between the 
diagrams in Figs. Qa) and Qb), we have 

S/oc = Sjj (6) 

if G is chosen such that 

G/oc — Gii (7) 

In standard fashioni, it can be seen from the Dyson equations for both the local problem 
and the lattice problem that, Eqs. ([],[?p are satisfied if the Weiss field is chosen as 

G 1 (iuj n )= iu n + /i - ^UotojlGijiiuJn) - G i0 (iuj n )G 0j (iuj n ) /G 00 (iuj n )} (8) 

u 

Consider now the higher-order correlation functions. In a skeleton expansion for any 
on-site correlation function of the lattice problem, every diagram also has the form of Fig. 
^|a) (with even(odd) number of fermion loops for any even(odd)-number-particle correlation 
functions). Likewise, the skeleton expansion diagrams for the corresponding local correlation 
function of the impurity problem have the structure of Fig. ^b). As a result Eqs. (§,0), or 
equivalently, Eqs. (£|,|8|), also guarantee that all the on-site higher-order correlation functions 
of the lattice problem can be calculated from the effective impurity problem. In particular, 

Xch,loc Xch,ii and Xs,loc ~Xs,ii- 
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Eqs. (§, H, f|) form the self-consistency equations. The effective impurity problem, 
defined by the action given in Eq. (|5]), can be equivalently written in terms of an impurity 
Hamiltonian of the following form, 



H imp = H kin + E d ^2 c\c a + Un cr n cl + tJ2( c lilka + H.i 
+Fj2-n c : {p g + pl q ) + + $- q ) 



q 



H kin = E k ril a ri k(T + ^ W q p\p q + ^ w q 4>\ ■ f q (9) 
ka q q 

where E d = —p and the parameters E k ,W q ,w q ,V,F, and g are given by 

iu n + p-t 2 J2l/(iuj n - E k )= G Q l {iu n ) 

k 

f 2 E wVtK) 2 - w 2 q ]= x£ Q (iu n ) 

q 

9 2 Y, W l/i( iu n) 2 - W 2 q }= X7,o(Wn) (10) 

<? 

where iu n and iv n are Matsubara frequencies for fermions and bosons respectively. 

Eq. (|9|) describes a single-impurity Anderson model coupled to two additional bosonic 
bands. The impurity corresponds to the c— orbital. r] kcr is the usual fermionic bath of the 
Anderson model, with a dispersion of E k . p q is a scalar-bosonic bath. <p q corresponds to a 
vector-bosonic bath. Note that the different components of the vector boson commute with 
each other: [0" , cjQ ■ ] = 8 a p8 qq i where a, (3 = x,y,z. The dispersions of the bosonic baths 
are W q and w q . 

As in the usual large D limit, having solved the local problem we can then calculate the 
lattice Green's function 

G(k,w) = — ^ (11) 

to + p — e k — S(o;) 

In the following, we establish that a parallel procedure can be carried out for the lattice 
correlation functions. 
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III. MOMENTUM-DEPENDENT CORRELATION FUNCTIONS 



We now specify the procedure to calculate the momentum-dependent correlation func- 
tions. For simplicity, we focus on the spin-spin correlation function only. The density-density 
correlation function has a similar form. 

A. Two-particle vertex functions 

We first establish the form of two-particle vertex functions. Consider the spin-spin cor- 
relation function, x s (q, uj), which can be written as 

X,(q, u) = J dexde 2 Xch{^ £2; q, w) (12) 
where €1,62 are illustrated in Fig. ||, which also specifies the Bethe-Salpeter equation, 

Xs,ij(e u e 2 ;uj) = x# (ei;o;)<y ei , ea + J dd 4imi( e i; w ) x 

Ii,mi,i 2 ,m2 

X^W2mim 2 ( e l> e 's UJ )Xs,l2jm 2 j( e > e 2~, (13) 

Here x ph represents the particle-hole bubble, with full fermion propagators. Ii 1 m 1 i 2 m 2 is the 
irreducible vertex function in the triplet particle-hole channel. It follows from the counting 
rules of Appendix |A| that only a limited number of contributions on the right hand side are 
of leading order. For the first term on the right hand side, only the i — j contribution is 
leading. For the second term, only terms with l± — rti\ — i and l 2 = m 2 contribute. Eq. 
(Pf) then leads to 

Xs{<h u)7i,e 2 = X ph (£i, ^)~ 1 ^ 1 ,, 2 - I{e u e 2 ; q, u) (14) 



where I(ei, e 2 ; q, u) = J2j eiq Rii hjij{ e ii e 2! Note that Eq. (|H|) is a matrix equation, with 
61,62 specifying matrix elements. 

We now need to calculate the irreducible vertex function I(ei, e 2 ; q, oS) in terms of the 
effective impurity problem. To do that, we first carry out a cumulant expansion for the spin- 
spin correlation function of the lattice model. A cumulant is introduced in a perturbative 
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expansion in terms of tij, Vij and Jij. It represents the corresponding local single-particle or 
two-particle Green's function calculated entirely in terms of the on-site ("atomic") part of 
the lattice Hamiltonian. For our purpose, it is more convenient to introduce an effective spin 
cumulant in analogy to the one-particle effective cumulant introduced by Metznei0. The 
effective spin cumulant is defined as all the diagrams for the on-site spin-spin correlation 
function which are irreducible in terms of cutting any J\ m line (where / and m are arbitrary). 
Loosely speaking, it is the bare spin cumulant plus all the local decorations. Denoting the 
effective spin cumulant as C s (ei, e 2 ; 00), the leading order spin-spin correlation function has 
the structure illustrated in Fig. |] and can be written as 

X;ij{ti, e 2] v) = C s (ei, e 2 ; w)<% + Y. de> de"C s {e 1 , e'; u) J u x,,y(e", £2; w) (15) 

1 J J 

which is equivalent to 

Xs(q,uj)-] £2 = C s (uj)-] e2 - J(q) (16) 

The leading in 1/D contributions to the effective spin cumulant contain diagrams of the 
form illustrated for the on-site self-energy in Fig. |2|. This again follows from the power 
counting rules of Appendix [A| and can be derived following a procedure parallel to that 
of Appendix [B]. As a result, all the loops involving intersite interactions are equal to the 
corresponding Weiss fields x^o an d Xcho- 

The above, in turn implies that the effective spin cumulant is equal to the local spin-spin 
correlation functions of the effective impurity problem, Eq. fl5|), which are not reducible in 
terms of cutting a single x7o li ne - We then have the following equation for the local spin-spin 
correlation function, 

Xs,ioc(e 1 ,e 2 ;uj) = C s (e u e 2 ; u) + J de' J de"C(e u e';uj) x^oM Xs(eV 2 ;^) (17) 
as illustrated in Fig. ||. Eq. (|T7|) leads to 

W-] e2 = xsM^)iL + x7,M (is) 



9 



which specifies how to calculate the effective spin cumulant from the Weiss field and the 
local spin-spin correlation function. 



Combining Eqs. (|TB , , we derive the following form for the momentum-dependent 



irreducible vertex function, 

J(ei, e 2 ; q, u)= x ph (e i; oj)- l 5 ei ,e, - X s,ioc^)lL ~ A^oM + A<l) (19) 
as well as the local irreducible vertex function, 

I Joc (ei, e 2 ; u) = x ph ^h w) _1 <J 6ll6a - Xs,ioc{u)~^ e2 ~ x7jo( w ) ( 20 ) 

B. Correlation functions 

We now proceed to determine the momentum-dependent correlation function. Again we 
focus on the spin-spin correlation function. 

Integrating both sides of Eq. fll5D over both ex and e 2 leads to the following, 

Xs,ij(u) = M s (u) S v + ]T M s {u)J ilXs ,ii{u) (21) 

i 

where the integrated spin cumulant M s {uj) = J dei j de 2 C s (ei, e 2 ; uS). Eq. (|2T| ) yields, 

M s (lo) 

Xs(Q,Uj) = ; ; ; : (22) 

The integrated spin cumulant M s (u) can be determined from integrating both sides of 
Eq. fllED over both e\ and e 2 , 



Xs,loc(u) = M s (uj) + M s (u) Xsfli^) Xb,Ioc( U ) ( 23 ) 

which yields, 

MM = XV ( \ ^ 

1 + Xs,o{u)Xs,loc{u) 



Inserting Eq. (El) into Eq. (Z2) leads to the final form for the momentum-dependent 



correlation function, 
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Eq. Q25D is one of the most important conclusions of this paper. It specifies how to 
calculate the momentum-dependent spin susceptibility from the local spin susceptibility, the 
Weiss field (x7,oi see Eq. (fD), and the exchange interaction. One check for the validity of 
this expression can be seen by re-writing it in the following form, 

x«(q,w) =Xs,ioc(u) + Xs,ioc(u)[J(<i) -x^oH]x»(<i» w ) ( 26 ) 

Using the expression for the Weiss field x7o> Eq. @, one finds that summing both sides 
of Eq. ( p6|) over q leads to Xs,ioc = Xs,ioc- An alternative derivation of Eq. fl25|) is given in 
Appendix y. 

IV. CONSERVING CRITERIA 

We see from the above that the single-particle self-energy is local, but the two-particle 
vertex function is momentum-dependent. We show in this section that, in spite of this, the 
approach is conserving. 

An approach is conserving if both the single-particle self-energy and two-particle irre- 
ducible vertex functions are variational derivatives with respect to the single-particle Green's 
functions of the Luttinger-Ward $— potential^. The latter is defined as the sum of all closed 
skeleton diagrams, and is a functional of the single-particle Green's function and interac- 
tion parameters. Applying this criterion to our case is somewhat subtle. The momentum- 
dependent part of the vertex function comes from diagrams for the Luttinger-Ward potential 
that are sub-leading in 1/D. These include the Hartree and Fock contributions. In fact, 
due to the usage of normal-ordered operators in the Hamiltonian the Hartree contributions 
to Luttinger-Ward potential and to the self-energy, given in Figs. |6]a) and |6|b) respectively, 
are identically zero. The corresponding Fock contributions, also given in Figs. ^|a) and 
|5|b) respectively, are sub-leading in 1/D. At the same time, both the Hartree and Fock 
contributions to the two-particle vertex functions, given in Fig. ||c), are of leading order. 



Because of this subtlety, we first use an alternate set of conserving criteria introduced 
in RefM Two conditions are sufficient for an approach to be conserving. The first states 
£2(1, 3, 1 + , 3 + ) = (^2(3, 1, 3 + , 1 + ), where G2 is the standard two-particle Green's function. 
[In this section we follow the short-hand notation of Ref.El. For instance, 1 labels (xi,n), 
and 1 + labels (x 1; ri + + ).] This condition is satisfied in our case. 

The second condition relates the single-particle self energy and the reducible two-particle 
vertex function. This condition contains two equations, one derived from the equation of 
motion, d T1 G(l, V) = [(7(1, 1'), H], and the other its adjoint. In what follows, we show that 
the first equation is satisfied; similar arguments lead to the validity of the second equation. 
In addition, for notational simplicity, we consider the case when only the intersite density- 
density interaction Vij is finite (the conclusion is unchanged if J^- is also present). In this 
case, the condition for the Hamiltonian given by Eq. ([]]) has the following form, 

£(1,10= -/ d3d5dQd8 v xljX _ [(7(1-3) 8(3 - 1') + (7(1, 5) (7(3,6) r(5,6,l',8) (7(8,3)] 
-JdESdS U [5(1- 1') n_ ff (l) + (7(1,5) (7(1,6) T(5, 6,1', 8) (7(8,1)] (27) 

where T is the reducible two-particle vertex function.0 

Consider first the case when xi is nearest neighbor to Xi'. The left-hand side is of order 
1/D and hence subleading. The right hand side contains the non-local component of the 
vertex function; one might then worry about whether it is still subleading. That it is so 
can be seen by enumerating the possible spatial locations of the integration variables and 
using the rules of Appendix [A|. In fact, the dominant contributions are all of order 1/D. We 
illustrate this with two examples. Consider first the case X3 = Xg = Xg = x 1; and x^ = x^, 
which contribute to the second term on the right hand side of Eq. fl27|). The non-local 
vertex function r(l', 1, 1', 1), which is of order l/\f~D, is accompanied by another non-local 
Green's function, (7(1, 1'), which is also of order 1/a/D. Consider next the contribution from 
the first term on the right hand side of Eq. (0), with X5 = xi' and X3 = Xg = x§ = x/, 
where x; is nearest-neighbor to Xj_. Again, the non-local vertex function r(l',Z, V,l), which 
is now of order 1/D, is accompanied by (7(1, 1') and v Xl)Xt , both of which are of order 1/yD; 
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summing over I gives rise to a factor D, leading to an overall 1/D contribution. We conclude 
that to leading order, Eq. fl2"T|) is indeed satisfied. Analogous arguments apply to arbitrary 

Xi ^ Xi/. 

Consider next xi = xi/. To the leading order Eq. ( p?D becomes, 

£(1,1) = ~Jd3 v XltX - G(l,l) G(3,3) r(l,3,l,3) G(3,3) 

-U [n_ CT (l) + G(l, 1) G(l, 1) r(l, 1, 1, 1) G(l, 1)] (28) 

[For notational simplicity, we have suppressed the time indices.] Here, the right hand side 
also contains the non-local components of the vertex function T. We can expand the right 
hand side in terms of the bare intersite interaction by expressing T in terms of the irreducible 



vertex function / and using Eqs. (19,2(j). This leads to diagrams which are in one to one 
correspondence with those for the on-site self-energy. Eq. fl28|) is then satisfied as well. 

What we have shown can be recast in terms of the derivability of both the self-energy 
and the irreducible vertex function. The only unusual feature is that, we need to keep not 
only the contributions to the Luttinger-Ward potential which are leading in 1/D, but also 
the Hartree-Fock-type terms which are formally subleading. The D — > oo limit has to be 
taken after variationally differentiating the Luttinger-Ward potential with respect to the 
Green's functions. The process of taking a variational derivative of a diagram with respect 
to the single-particle Green's function can change the order in 1/D of its contribution. 

V. DISCUSSIONS 
A. Multi-band models and ordered states 

The generalization of our approach to multi-band cases is straightforward. For a two- 
band extended Hubbard model, the corresponding mean field equations have already been 
written down in Refs.@@. 

Another important two-band model is the Kondo lattice model, 
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H = tijCiaCjv + Yl J K$i ' + J H^i ' $3 ( 29 ) 

<ij>,a i <ij> 

where Si denotes an impurity spin at site i, and s Ci represents the spin of conduction (c— ) 
electrons at site i. Taking the large D limit, again with = to/y/D and Jy = Jq/\/T), 
results in the following effective impurity action, 

S MF = S top + f drJ K S-s c 
Jo 

- fdr r^Ect(r)G^-O^(r0 + ^(r)-^(r -0^(70] (30) 
Jo Jo „ 

where S top describes the Berry-phase of the impurity spin. The Weiss fields Gq 1 and x^o 
are determined by the self-consistency equations as given in Eqs. (HU). The effective action 
can equivalently be written in terms of the following impurity problem, 

Himp = E kvLvk* + Y W A ■ <S? ~ V- c l c ° + tJ2( c lVk* + H.c.) 

fccr <? O" fccr 

+J K S-s c + gJ2S-($ q + $l q ) (31) 
q 

where Ek, t, w q and g may be determined from the Weiss fields Gq 1 and Xs,o as specified by 
Eq. 

Finally, we can also extend the approach to a state with long-range commensurate spatial 
ordering. This requires taking the normal-ordering, as specified in Eq. (ffl), with a site- 
dependent average charge or spin appropriate for the ordered state. The dynamical mean 
field equations, Eqs. (|||| §) still apply. 



B. Incommensurate susceptibilities 



The form of the momentum-dependent correlation function given by Eq. (25) applies 
to generic q. The momentum-dependence is entirely given by that of the corresponding 
intersite interaction. It does not depend on the single-particle dispersion. 

The situation is different from the standard large D limit, where the momentum depen- 
dence of the correlation functions is given entirely by the single-particle dispersion. The 
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latter is possible for lattices with unbounded bare density of states - such as hypercubic 
lattice - which contains special q such that e q = J2ij e iq,Rii tij is of order \I~D. 

Formally, the extended DMFT described here can only be defined for lattices with a 
bounded bare density of states (such as Bethe lattice). It becomes ill-defined for lattices 
with unbounded density of states. When e q is of order Vd, so is J(q); through the form of 
the susceptibility, Eq. fl2~5[) , the system would then become unstable (when J(q) is positive). 

One way to approximately incorporate the incommensurate fluctuations induced by 
Fermi-surface features is through the (exact) Bethe-Salpeter equation 

Xs (q, u)-] e2 = X p/l (ei; q, - I(e 1 , e 2 ; q, u) (32) 

where x ph { e u <L u ) is the usual particle-hole susceptibility bubble calculated in terms of 
the full single-particle propagators. One uses the the self-energy and the irreducible vertex 
function, Eq. ([19]), using the extended DMFT. At the same time, one uses the momentum- 
dependence of the intersite interactions in a given system at finite dimensions. This proce- 
dure is of course no longer systematic. 

C. Comparison with other approaches 

A direct 1/D expansion has been introduced by Schiller and IngersentB The expansion 
is carried out for the Luttinger-Ward $— potential^!, i.e., all the closed skeleton diagrams. 
The leading (1/D°) order contributions involve, as usual, only a single site. The next-to- 
leading order (1/D) contributions come from diagrams involving two sites. A dynamical 
mean field description to this level requires two effective actions describing a single impurity 
and a two-impurity cluster, each coupled to its respective self-consistent medium. Similarly, 
expanding to order 1/D n requires solving simultaneously effective problems involving one- 
site, two-site, and up to n + 1 site clusters, each embedded in its own self-consistent medium. 

The extended DMFT described here can be thought of as a conserving resummation of 
contributions to all orders in 1/D. This is illustrated in Fig. (|7|). A detailed comparison 
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between the two approaches should in principle be meaningful when spatial correlations are 
short-ranged. 

An alternative approach has recently been introduced by Hettler et alM In this ap- 
proach, the Brillouin zone is divided into N c regions. For each frequency, one introduces 
one Weiss field for each momentum region. The resulting self-consistent problem describes 
an N c — site cluster embedded in these N c self-consistent Weiss fields. For iV c = 1, it reduces 
to the standard large D DMFT. The precise relationship between their approach and the 
extended DMFT described here is unclear. 

Finally, our approach also has some similarities with the dynamical mean field theory of 
random spin systems. Bray and Moore@ considered the quantum Sherrington-Kirkpatrick 
model, in which the exchange coupling is infinite-ranged and has a Gaussian distribution of 
mean zero and variance scaled to J 2 /N where N is the size of the system. Carrying out 
disorder averaging using replicas, as usual, leads to a problem with four-spin interactions; 
each spin now carries a replica index. Since the exchange couplings associated with different 
bonds are uncorrelated, the four spin interaction has the form of two spins at different times 
t and t', from any site, interact with two other spins, also at r and r', at every other 
site. Taking the N — > oo limit then leads to a single-site problem with retarded spin-spin 
interactions. In the paramagnetic phase, the effective impurity problem has a similar form 
as the spin part of Eq. (|5]), with a self-consistency equation which is also similar - though 
not identical - to Eq. (^). Generalizing^B to the case when conduction electrons are 
also present gives rise to mean field equations similar to those described here. The form of 
correlation functions is of course very different. 

VI. CONCLUSIONS 

To summarize, we have given a perturbative derivation of an extended dynamical mean 
field theory. This approach goes beyond the standard D = oo DMFT by incorporating the 
quantum fluctuations associated with intersite interactions. 
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The self-consistent impurity problem has the form of an Anderson impurity model with 
additional bosonic baths. This is a novel kind of impurity problem, and is of interest in 
its own right. The role of a scalar boson bath is to introduce additional screening, thereby 
enhancing orthogonality. The precise consequences of the enhanced screening depends on the 
form of the spectral function of the bosons. Ref.B analyzed an impurity problem containing 
both a scalar boson bath with a spectral function of u A and a conduction electron band with 
a regular density of states at its Fermi energy. The correlation functions have mean-field 
exponents. The effect of an anisotropic vector-boson bath is similar to the scalar case. The 
effect of an isotropic vector-boson bath is more complex. We and independently Sengupta 
have carried outi'0 a renormalization group analysis of spins coupled both to a regular 
conduction electron band {Jr) an d to a vector-boson bath (g) with a spectral function 
of u A . In the sub-ohmic case (A < 1), there exists a critical point separating a phase 



characterized by the fixed point at J K = oo,g* = and another by J K = 0,g* ~ yl — A. 
The critical exponents are anomalous^. In the absence of conduction electrons, the impurity 
problem corresponds to a spin coupled to vector bosons alone. Sachdev and Ye§l solved such 
a model in the large N limit (see also Ref.il). Finally, we note that quantum Monte Carlo 
methods have recently been developed for this type of impurity problems^. 

The momentum-dependent correlation functions in the extended DMFT have the general 
form given by Eq. (f25|). While different in details, it is similar to the phenomenological 
expression recently introduced in Refs. in the context of the dynamical spin susceptibility 
in a heavy fermion metal {CeCu^ x Au x ) close to a zero-temperature phase transition. The 
experiments§lizl at the critical concentration are suggestive of a dynamical spin susceptibility 
[J(q) + /(u;,T)] _1 over the entire Brillouin zone, where the frequency and temperature 
dependence in f(u,T) show an anomalous exponent. Whether a quantum critical point 
with an anomalous exponent in the dynamical susceptibility can emerge in the extended 
DMFT described here is an exciting open question. 

Finally, the extended DMFT described here should be applicable to quantitatively study 
the short, but finite, ranged dynamical fluctuations in paramagnetic phases (whether or not 
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the ground state is ordered). Consider, for example, the spin fluctuations in the paramag- 
netic phase of a Mott insulator. In the standard large D limit, these fluctuations reflect 
simply isolated local moments: They contain no damping and are featureless in momentum 
space. In the extended DMFT, the self-consistent problem introduces mode coupling and 
hence damping. In addition, through Eq. Q25D the dynamical spin susceptibility is expected 
to be peaked at antiferromagnetic wavevectors. This approach is particularly useful at tem- 
peratures where the correlation length is short; here approaches taking into account only 
long-wave-length fluctuations would break down. Quantitative calculations using the ex- 
tended DMFT will perhaps allow a detailed understanding of the neutron scattering results 
in both undoped and doped Mott insulators, such as V2O3S where the exchange interactions 
are not particularly large. 
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APPENDIX A: POWER COUNTING RULES 



The order in 1/D of the correlation functions and vertex functions can be determined 
by analyzing the diagrams in real space. 

Consider first the single particle Green's function and vertex functions Iijij and Fyy. 
In any real space Feynman diagrams for these quantities, it takes at least ||i — j\ \ number 
of fermion propagation or intersite interaction steps. Therefore, Gij ~ i^y ~ r^y ~ 
0(i)IIMI/2. 

Consider next the vertex functions involving three independent sites, I(i,j,i,l) and 
T(i,j,i,l). Every diagram involves at least three independent paths of fermion propa- 
gators and/or interaction lines. Therefore, ~ T{i,j,i,l) ~ o 

(l/£))^8/2 where 

M 3 = min (||z - j|| + + lb'-i||,||*-z|| + 

Finally, consider the vertex functions involving four independent sites, I(i,j,l,m) and 
T(i,j,l,m). Every diagram involves at least four independent paths of fermion propagators 
and/or interaction lines. As a result, I(i,j,l,m) ~ T(i, j,l,m) ~ o(-^) M4/ ' 2 where M 4 = 
min(Xi,X 2 ,X3,X 4 ). Here, Xx = \\i-j\\ + \\i-m\\ + \\l-m\\, X 2 = \\j-i\\ + \\j-l\\ + \\j-m\\, 
X3 = \\l — i\ \ + \\l — j\\ + \\l — m\\ and X4 = ||m — z|| + ||m — j\\ + ||m — 

APPENDIX B: ON-SITE SELF-ENERGY 

Consider an arbitrary skeleton expansion diagram for the on-site self-energy for site 0, 
Eoo(^)- From the counting rules of Appendix |A], only local fermion propagators appear. 
The only non-local terms involve intersite interactions. These intersite interaction terms 
can be grouped into separate loops, each starting at site and returning to site as is 
illustrated in Fig. 0a). Note that, no loop can return to site more than once. For instance, 
a contribution given in Fig. [8] is subleading. 

The analytic expression for each loop can be determined as follows. The beginning and 
ending interaction lines give a product JoiJ m o, where / and m are arbitrary sites nearest- 
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neighbor to site 0. The solid square then represents a correlation function involving Si and 
S m . Given that site is excluded from anywhere in the solid square, this correlation function 
has to be evaluated in terms of H^, defined as the original Hamiltonian Eq. ([]]) with site 
excluded. As a result, this loop can be written as 



,(o) 

Im 

,(°) 



^JolJmOXlm (Bl) 



where xim =< T t Si(t) ■ S m (r') > H (o), which can be determined as follows. For Xim =< 
T t Si(t) ■ S m (r') > H , the cumulant expansion given in Eq. flT5|) implies \im = XiiXimXmm, 
where x'lm = T, pa ths JihXhhJhhXhh-Xi n i n Ji n m and [/, h, l 2 , l n , m] labels a non-self- 
retracing path from site / to site m. This in turn implies Xim = Xim ~ XuXioXooXomXmm- 
Therefore, 

Xim = Xim - Xl0X0m/X00 (B2) 

leading to the expression for the spin Weiss field, Xs,o> given in Eq. (f|). 

Similarly, an interaction chain generated by t>y corresponds to the charge Weiss field, 

Xrf£o» s iven in Ec i- (D- 

APPENDIX C: ALTERNATIVE DERIVATION OF THE 
MOMENTUM-DEPENDENT SUSCEPTIBILITY 

In this section we present an alternative derivation for the momentum-dependent sus- 
ceptibility x(q, w), Eq. (H). 

We re-write Eq. @ in momentum-space, 

xo(^)- 1 = E AoMo, <*>) - E J(q)x(q, ")?/xioM (ci) 

In addition, we use Eq. (E^), i.e., 

»^'- M(,)-'- JW < C2 > 



We now substitute Eq. (|C|) into Eq. (|GT[) . By using 
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and 

^ MM J ' q -J( q ) ^ M( " ) - 1( - 1 + ^ )) (C4) 

we obtain XoM -1 = M(uy l - xiod^y 1 , i.e., 

mh-^x^M^ + xoH" 1 (C5) 

Inserting Eq. flC5|) into Eq. Q then leads to Eq. (H). 

Note that, Eq. ([25|) reduces to the correct result for generic q in the standard large D 
limit where J(q) = and Xo 1 = 0- 
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FIGURES 




FIG. 1. A typical intersite self-energy diagram. A full line represents a single-particle Green's 
function, and a dashed line corresponds to an intersite interaction. This diagram is sub-leading. 




FIG. 2. a) An on-site self-energy diagram of the lattice problem, b) The corresponding local 
self-energy diagram of the effective impurity problem. A shaded line corresponds to a Weiss field. 
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FIG. 3. The Bethe-Salpeter equation. Here a bubble with a shaded insertion represents the 
susceptibility. The shaded square describes an irreducible vertex function. 
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FIG. 4. Susceptibility expansion in terms of cumulants. The bubble with a shaded insertion 
represents the susceptibility, while the bubble with a solid insertion corresponds to an effective 
two-particle cumulant. 
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FIG. 5. Local susceptibility expansion in terms of cumulants and the Weiss field. 
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FIG. 6. Hartree and Fock contributions to a) the Luttinger-Ward potential; b) the self-energy; 
and c) the vertex function. 
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FIG. 7. a) single-site, two-site and three-site diagrams for the Luttinger-Ward potential in a 
direct expansion in 1/D; b) The corresponding diagrams in the extended DMFT. 
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FIG. 8. A diagram for the on-site self-energy. It connects the selected site to the rest of the 
lattice by three interaction lines, making it subleading. 
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